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Abstract. The behaviour of periodic points of discrete Euler top is studied. We derive 
| invariant varieties of periodic points explicitly. When the top is axially symmetric they 

are specified by some particular values of the angular velocity along the axis of symmetry, 
£C} , different for each period. 
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Ph; 1 Introduction 

J: 

The Kowalevski workshop on mathematical methods of regular dynamics was organized by 
Professor Vadim Kuznetsov in April 2000 at the University of Leeds [I] . In his introductory talk 
about the Kowalevski top, Professor Kuznetzov [2] had shown his strong interest on the subject 
. £h \ and motivated the authors to work on classical tops. 

^ ' In our recent paper [3] we have studied the behaviour of periodic points of a rational map 

and found that they form a variety for each period specified by invariants of the map if the map 
is integrable, while they form a set of isolated points dependent on the invariants otherwise. 
It is apparent that an application of our theorem to the problems of a classical top is quite 
interesting and will be fruitful. We investigate the discrete Euler top, in this article, to see how 
the invariant varieties of periodic points look like in this particular example. In conclusion we 
will show that there is no periodic points of period 2 and 4 if the top is not axially symmetric. 
In the case of period 3 we derive explicitly an algebraic variety of dimension two as an invariant 
variety of periodic points. When the top is axially symmetric, the angular velocity of a periodic 
map along the symmetry axis is quantized to some special values determined by the period 
and the shape of the top. The other components of the angular velocity are free, thus form an 
invariant variety of periodic points separately for each period. 

To start with let us briefly review our theorem of [3]. We consider a rational map on C d , 
where C = {C,oo}, 

x = (xi,x 2j ...,s d ) - X = (X 1 ,X 2 ,...,X d )=:X( 1 \ (1) 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics". 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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We are interested in the behaviour of the sequence: x — ► — ► X^ 2 ) In particular we 

pay attention to the behaviour of periodic points of rational maps. If the map is nonintegrable 
we shall find a set of isolated points with fractal structure as a higher dimensional counterpart 
of the Julia set. Our question in [3] was what object appears when the map is integrable. 

We assume that the map has p (> 0) invariants Hi(x), H2(x), . . . H p (x). If hi, h,2, ■ ■ ■ , h p 
are the values of the invariants given by the initial values of the map, the orbit of the map 
is constrained on the d — p dimensional variety determined by the conditions Hj(x) = hi, 
i = 1, 2, . . . ,p, which we denote by V(h). Now let us consider the periodicity conditions X( n ) = x 
of period n. We can eliminate p variables out of x and only d — p periodicity conditions remain. 
If they are independent, we obtain a set of isolated periodic points on V(h) in general. 

It may happen, however, that some of the d—p conditions impose some relations on hi, h%, ■ ■ ■ , 
hp instead of fixing all of the d—p variables. If m is the number of the conditions which determine 
the values of the variables, d — p — m variables are left free and the periodic points of period n 
form a subvariety of dimension d — p — m on V(h), instead of a set of isolated points. In this 
case we say that the periodicity conditions are 'correlated'. If / is the number of the periodicity 
conditions which relate the invariants, p — I invariants remain undetermined. This means that 
the periodic points of period n form a d — I — m dimensional subvariety in C d . We have proven 
in [3] the following lemma: 

Lemma 1 ([3]). A set of correlated periodicity conditions satisfying mm{p, d — p} > I + m 
and a set of uncorrelated periodicity conditions of a different period do not exist in one map 
simultaneously. 

When m = the periodicity conditions determine none of the variables but impose I relations 
among the invariants. In this particular case all points of V(h) are the points of period n, 
while the variety V(h) itself is constrained by the relations among the invariants. We call the 
periodicity conditions are 'fully correlated' in this case. The periodic points form a subvariety 
of dimension d — I in the space C d , which we call 'an invariant variety of periodic points'. Every 
point of this variety can be an initial point of the n period map, whose orbit stays on it. Since 
the condition minjp, d — p} > I + m is automatically satisfied, our theorem follows to the 
Lemma Q] immediately: 

Theorem 1 Q3J). If there is an invariant variety of periodic points of some period, there is no 
set of isolated periodic points of other period in the map. 

The Theorem [1] doesn't tell us directly whether the map is integrable or nonintegrable. There 
is, however, some evidence to believe that the periodic points of a nonintegrable map, if they 
exist, form a fractal set of isolated points. Therefore it is reasonable to adopt the following 
proposition as our working hypothesis: 

If a map is nonintegrable, there is a set of uncorrelated periodicity conditions of some period. 

The Julia set, which is the source of chaotic orbits, is a subset of the closure of all isolated 
periodic points. We emphasize that our hypothesis does not require that all of the periodicity 
conditions of a nonintegrable map are uncorrelated but requires only one at least. On the other 
hand our theorem shows that the existence of an invariant variety of periodic points excludes a set 
of isolated periodic points in the map and vice versa. This means that if a set of periodic points 
of some period forms an invariant variety there is no Julia set, thus suggesting the following 
statement: 

Conjecture ([3]). If there is an invariant variety of periodic points of some period, the map is 
integrable. 

Note that this does not exclude possibilities that some integrable maps do not have an 
invariant variety of periodic points. For example an integrable map with no invariant does 
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not have an invariant variety. On the other hand there are d dimensional maps which have 
d — 1 invariants but reduce to the logistic map after elimination of d — 1 variables by using 
the invariants. Therefore the situation is quite different from the continuous time Hamiltonian 
flow, whose integrability is guaranteed by the Liouville theorem if there are sufficient number 
of invariants. 

In order to support our conjecture we have studied in [3] various maps, such as the QRT 
maps, the Lotka-Volterra maps and the Toda maps, and found that all periodic points form 
invariant varieties of periodic points if the map is integrable and periodic points exist, while 
no such property has been found otherwise. We also studied the g-Painleve maps which are 
integrable but not volume preserving in general. We found invariant varieties only when the 
parameters are restricted so that the maps have sufficient number of invariants. 

We introduce the discrete Euler top in Section 2. To find an invariant variety of periodic 
points for the Euler top the general scheme developed in [3] is applied in Sections 3 and 4. In 
the final section we study explicitly the nature of the invariant surfaces of an axially symmetric 
Euler top. 



2 Discrete Euler top 

When the time is continuous, the equation of motion for the Euler top is given by 

T duJi . , du 2 ,j T s T duj 3 

dt dt dt 

where (ui, u 2 , u 3 ) are the angular velocity in the body fixed frame and (ii, I21H) are the cor- 
responding moments of inertia. The system has two invariants, ^(I\ujf + l2^\ + ^3^1) an d 
Iiiuf + I 2 w\ + corresponding to the total kinetic energy and the square of angular mo- 

mentum, hence is integrable. 

A discretization of the Euler equation, which preserves integrability, was first obtained by 
Bobenko et al. in [3], and then discussed by other authors [5J El [8] . We adopt here an explicit 
version of the discretization proposed by Hirota et al. [3 [6]. After the discretization we write 
the angular velocity as (xi,X2,X3) instead of (ui\, ui 2 , w 3 ) and consider the map (xi,X2,x 3 ) — ► 
(X 1 ,X 2 ,X 3 ) defined by 

h{X x - Xl ) = ^{h~ h)(X 2 xz + x 2 X 3 ), 

h(X 2 - x 2 ) = \(h~ h)(X 3 xi + x 3 Xx), (2) 
I 3 (X 3 - x 3 ) = -{h - h){X lX 2 + xxX 2 ). 

The continuous limit corresponds to 5 — ► 0. 

Solving the equation ([2]) for (X±, X2, X 3 ) we find 



X-i 



x\(l — a,20L 3 x\ + a 3 a\x\ + a±a2X^) + 2a±X2X 3 

1 — 2010203X1^2X3 — 02Q3xf — 0301X3 — Oi02x|' 

x 2 (l + 0203X1 - 0301X3 + 0102X3) + 2o 2 x 3 xi 

X<2 ~ \ ^ 2 2 2> V>) 

1 — 2010203X1X2X3 — 0203XJ — 0301X3 — 0102X3 

X3(l + 0203X1 + 0301X2 — 0102X3) + 203XlX 2 

A3 - 



1 — 2010203X1X2X3 — 02Q3xf — OL 3 OL\x\ — 0102X3' 

where we used the notations 

I2 — I3 I3 — h h — I 2 
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We notice that, when the top is axially symmetric, the map ([3j) is nothing but a two dimen- 
sional linear transformation. For example if we assume I2 = I3, the map becomes 



X2 = X2 cos Q, + X3 sin 0, 
X3 = X3 cos f2 — X2 sin £1, 



(4) 



where 



cos 



4g - (h - hfxj 
4/ 2 + (/ 2 -/i) 2 * 2 ' 



Therefore we discuss the axially symmetric top separately from the generic case in Section 5. 
The invariants of the map ([3]) are 



I\x\ + hxl + hx 



1 — 0L2Ci-i,x\ 



Ho 



l\x\ + + I 



2 T 2 
3 X 3 



1 — Ol2Ci-iX\ 



(5) 



as it will be checked by a direct substitution of ([3j). They coincide with the invariants of the 
continuous case in the limit 5 — > 0. We fix the value of (5 at 1 hereafter, since it is irrelevant in 
the following discussions. 

If we denote by x one of the three variables (x\, X2, £3), the elimination of other two variables 
from the map (|2|) yields 



S(X,x;q) = 

with 

S(X, x; q) := aX 2 x 2 + bXx(X + x) + c(X - x) 2 + dXx + e(X + x) + f. 
The parameters q = (a, b, c, d, e, f ) are given by 



(6) 



(7) 



/ 


«1 








bi 








ci 








di 








ei 






V 


/1 








r/ 2 








&2 








C2 








d 2 














V 


/2 







\ 



/ «3 \ 

63 
C3 

d 3 
e3 
V /a / 





(A - A 2 + A3) 2 

4(a 2 - A 2 (A + A 3 ) + A 3 (A - A 2 ) 


V (4/0203)^2^3 

-4o 3 ai A (A - A 2 ) 


(A) -M- A3) 2 

4(A 2 -A (A2+A3)-A 3 (Ao-A2) 



V (4/0300^3^1 

/ -4oia 2 Ao(A + A 3 ) 


(A + A 2 + A3) 2 

4(A 2 + A (A 2 + A 3 ) + A 2 (A + A 3 ) 


\ (4/oio 2 )AiA2 



(8) 



\ 
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corresponding, respectively, to x = x±,x 2 ,xs. Here we introduced the notations 

A) = 4/i/ 2 / 3 , Ax = (I 2 - h){hHx - H 2 ), A 2 = (J 3 - h)(I 2 H x - H 2 ), 
A 3 = (h - h)(hHx - H 2 ), (Ax + A 2 + A 3 = 0) . 

One might wonder that the map x — ► X defined by ([6]) does not determine an image of 
the map uniquely. Since the function S(X, x; q) of ([7]) is symmetric under the exchange of 
the variables x and X, we see that the two solutions of © correspond to the forward and the 
backward maps, which we denote and X^. If we apply the map to X^ we should get x 
and X( 2 \ In this way we shall obtain a chain of images of the map into two directions: 

. . . ^_ X (~ 2 ) <— X^ ^x^ X® — X( 2 ) 

The number of branches of the map does not increase, but remains two, at every step of the 
map. Therefore a map of the form ([6]) is well defined in general. 



3 Iteration of the map 



We studied in [3] the map defined by ([6]) and called it a 'biquadratic map'. An iteration of the 
map yields the biquadratic map again but with new parameters q( 2 ); 

(2) 



a 

6 (2) 



(ae - cb) 2 - (ad - 2ac - b 2 )(be -cd + 2c 2 ), 
(ae - cb)(2af -be + cd- Ac 2 ) - (ad - 2ac - b 2 )(bf - ce), 
(af-c 2 ) 2 -(ae-bc)(bf-ce), 

4(a/ - c 2 ) 2 - 2(ae - 6c) (6/ - ce) - (6e - cd + 2c 2 ) 2 (9) 
- (ad - lac - b 2 )(df - 2c/ - e 2 ), 
e ( 2 ) : = (fb - ce)(2a/ - fee + cd - 4c 2 ) - (/d - 2/c - e 2 )(ea - c6), 
/( 2 ) : = (fb - ce) 2 - (fd - 2/c - e 2 )(fee - cd + 2c 2 ). 

If we repeat the map further we obtain a series of biquadratic maps whose parameters can 
be determined iteratively from the previous ones as follows: 



>+i) 



1 



^((a A c) 2 - (a A b) n (b A c) n ), 



1 



<7 



(n-1) 



pr((a A fe)„(fe A c) n - (a A c) 2 ) + (a A c) n ,((a A e) n + 2(fo A c)„) 



>+i) 



- 2 (( a A fe )n( 6 A e )i ~ (a A fe) n (c A d) n + (a A d) n (6 A c) n ) J , 

_L_(( ce M _ 6/M)( ae (n) - 6c (n) ) + (cfe (n) - ea (n) )(/fe (n) - ec (n) ) 
+ (a/W - cc (n) ) 2 + (/a (n) - cc (n) ) 2 ), 



(10) 



= _L_ ( _ / (n-l) a (n+ 1 ) _ a (n-l) ; (n+l) _ ^(n-iyn+l) _ ^(n-lj^n+l) + (a/J) 2 

+ ( c Ad)l - (a A b) n (e A f) n - (b A c) n (c A e) n + (a A d) n (d A f) n + 2(b A e) n (a A f) n 
- (3(c A e) n - (b A f) n - (d A e) n ) (3(b A c) n - (a A e) n - (b A d) n ) 

+ 2((a A d) n - (a A c) n ) ((c A f) n - (d A f) n ) + 2((b A c) n + (a A e) n ) ((b A f) n + (c A e) n fj 
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e(n+1) = 7(^1) ( 7^=i) ((/Ae)n(e A c)„ - (Ac)*) + (/ A c)„((/ A 6)n + 2(e A c) n ) 
1 



p(n+l) 



((/Ae)n(e A &)n - (fAe) n (c A d) n + (/ A d)n (eAc)n) 



j (n _i) ((/ac)» - (/ A e)n(e A c)„), 



where we used the notation (p A g') n = gg'^ — g'g^ ■ 

Despite the complicated expression of the relation (fTUj) , we observe a special dependence on 
the nth parameters q( n ). Besides c^ n+l \ the dependence of the (ra + l)th parameters on the rath 
ones is always in the form (g^g^n = gg ,{xL) - g'g^ 1 ' '• They all vanish simultaneously when the 
parameters q( n ) are 'fully correlated', that is, if there exists a function 7^ n+1 ^(q) such that 

q (n) =eq + 7 (n+1) (q)q (n) , (11) 
where e is an arbitrary constant. In fact we obtain, after some manipulation, 

(a A b) 2 = (af -eb- 3c 2 + cd)(2a 2 e - abd + b 3 ), 
(a A c)2 = (af — eb — 3c 2 + cd)(a 2 f + ac 2 — acd + b 2 c), 
(fr A c)2 = {af — eb — 3c 2 + cd)(2ace — abf — be 2 ), 



(e A /) 2 = (af -eb- 3c 2 + cd)(e# - e 3 - 26/ 2 ), 
from which we find 7^(q) : 

7 (3) ( q ) = a f -be- 3c 2 + cd. (12) 
The formula (|1Q|) enables us to find a series of (q) iteratively, as follows: 

7 (4) (q) = 2ac/ - adf + b 2 f + ae 2 - 2c 3 + c 2 d - 2bce, (13) 

7 ( 5 )(q) = a 3 / 3 +(- cf 2 d + 2cfe 2 + fde 2 - 3ebf 2 - e 4 - c 2 / 2 )a 2 

+ ( - 13c 4 / + 18c 3 /d + de 3 b + 2cf 2 b 2 + 7dc 2 e 2 - ce 2 d 2 - 2ce 3 b 

+ 2c 2 feb - 7fd 2 c 2 - Uc 3 e 2 + cd 3 f + fb 2 e 2 + f 2 db 2 - ebd 2 f^j a (14) 

- cd 2 b 2 f - b 3 e 3 - Ac 3 deb + cdb 2 e 2 + 13ec 4 6 - f 2 b 4 + 7fb 2 c 2 d 
+ c 4 d 2 - 5c 5 d + 5c 6 - 2fb 3 ec - e 2 c 2 b 2 + eb 3 df - Ufb 2 c 3 , 

and so on. 

When (fTTj) holds, the equation S(Q, x; q n +i) = can be written as 

c («+l)(Q - x) 2 + ( 7 (n+1) (q)) 2 ^n+i(Q,x) =0, n = 2,3,4, ... . (15) 

Here 

K n+1 (Q,x) = a^Q 2 x 2 + U n+l \Q + x)Qx + d {n+ ^Qx + e {n+l \Q + x) + f n+1 \ 

and a^ n+1 \ for instance, is obtained from a( n+1 ) simply replacing (g/\g') n by (<7 A g')n- If Q is 
a point of period n + 1, the first term of (|15p vanishes. Hence the periodicity condition requires 
for the second term to vanish. This is certainly satisfied for arbitrary x if 7^ n+1 ^(q) = 0, namely 
when the periodicity conditions for the parameters q( n ) are fully correlated. The other possible 
solutions obtained by solving K n+ \(x,x) = do not correspond to the points of period n + 1, 
but represent the fixed points or the points of periods which divide n + 1. 
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4 Invariant varieties of periodic points 
for the discrete Euler top 

We are ready to study the periodicity conditions for the discrete Euler top. Throughout this 
section we will not consider axially symmetric cases, which we discuss in the next section. The 

(n) 

direct calculation of the periodicity conditions X- = Xj is not easy to carry out by a small 
computer. Therefore we use the method we developed in the previous section. Before starting, 
however, let us first search the fixed points of the map. If we remember that the variables 
(xi, x 2 , x 3 ) are the discrete analog of the angular velocity (ui, u 2 , u 3 ), a fixed point of the map 
corresponds to the motion of the top which does not change the angular velocity in all directions 
of the body fixed frame. 

Needless to say the fixed points are nothing to do with the invariants of the map. To find 
them we go back to the map ([2} and see immediately that they are 

fixed points : {x | X2 = x 3 = U x 3 = x\ = U x\ = X2 = 0}. 

This result shows that a fixed point is realized as a steady rotation along one of the three axes. 
The value of the angular velocity along the direction is arbitrary, while the angular velocities 
are zero along the other two directions. 

The method we developed in the previous section enables us to get information of the period- 
icity conditions of period greater than three. The case of period 2 must be considered separately. 
From the general expression (|9|) the parameters in S(X^ , x; q 2 ) are not proportional to a com- 
mon factor. But if we substitute ([8|) into (|9|), we find that they have the following form 

a< 2 > = ( 7 (2) (q)) 2 « (2) , & (2) = 0, = ( 7 {2) (q)) 2 ^ 2 ), 

e (2) =0, /( 2 ) = ( 7 {2) (q)) 2 / (2) 



with 

7 i 2) (q) = A -A 2 + A 3 , 7 f ) (q) = A - A 2 -A 3 , 7 J 2) (q) = Aq + A 2 + A 3 , (16) 



corresponding to x = X\, x 2 


, x 3 , 


respectively. For the higher periods we can apply the formu- 


lae {32]) and 


(fl~3l) to obtain 








(3) 

7i = 


(A 2 + 2A (A 2 - 


A 3 ) 


-3A 2 )((A 1 + A ) 2 + AA A 3 ), 




(3) 
7 2 = 


(A 2 1 + 2A (A 2 - 


A 3 ) 


-3^ 2 )((^ 1 + A ) 2 + 4A 3 ^i), 


(17) 


(3) 
7 3 = 


(Aj + 2A (A 2 - 


A 3 ) 


-3^ 2 )((^ 1 -A ) 2 + 4A 1 ^ 2 ), 




(4) 

7i = 


2(A 1 -A )(A - 


-A 2 


- ,4 3 ) ((^i - A)) 4 - 8A 2 (A + A 3 )(A\ + A 2 )) , 




(4) 
7 2 = 


2(A -A 1 )(A - 


-A 2 


+ A 3 ) ((A! + A ) 4 + 16A AxA 3 (A - A 2 )) , 


(18) 


(4) 
7 3 = 


2(A + A 1 )(A - 


-A 2 


+ A 3 ) ((Ax - A ) 4 + 16A AxA 2 (A + A 3 )) , 





From the expressions (116|) . fjlTf) . (I18p it is clear that the periodicity conditions do not deter- 
mine points but impose relations among the invariants of the map. This owes to the fact that 
the initial parameters q are dependent on the invariants of the map, as we see in 0. After 
iteration of the map n times the new parameters q( n ) are also dependent on the invariants. 
Therefore the periodicity condition 7 ^ n ^(q) = imposes relations among the invariants. 
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The periodicity conditions of period n are satisfied when 7} = 0, 72 = 0, 73 = are 
satisfied simultaneously. By an inspection of (|16p . (|17p . (|18p we notice that the conditions are 
satisfied by the single condition 

Aj + 2A (A 2 - A 3 ) - 3A 2 = (19) 

in the case of n = 3. All other cases impose further relations among the invariants. 

In order to find where the periodic points are in the (xi, X2, £3) space, we simply substitute 
the formulae © into the conditions 7^ (q) = 0. In the case of (|19|) we obtain 

= {x| (1 + 6+6 + 6) 2 - 4(1 + 66 + 66 + 66) = 0} , (20) 

in terms of the new variables 

= ^ ^ ^ = ^ x 2l 6 = ^ x 3 . 

The set of points satisfying (|20|) form a variety of periodic points of period 3 in the space of 
(x±, X2, X3). Every point on this variety is a point of period 3. We called this type of variety 
'an invariant variety of periodic points', because it is determined uniquely by the invariants of 
the map alone. The dimension of the variety is two, which is the number of the invariants. In 
the case of (120p the invariant variety is an algebraic variety of degree 4, symmetric in the three 
variables X\, X2, £3- 

Now let us pause a while. The discrete Euler top (J2|) has been known being satisfied by elliptic 
functions as special solutions. The map generates an elliptic curve. This curve, however, is not 
the invariant variety of periodic points in our consideration, since the map is not controlled, in 
general, by the periodicity conditions of some fixed period. In fact the invariant variety of (|20p 
is not a curve but a surface. Once an initial point is chosen on the surface, the orbit stays on 
it before it returns to the initial point. The invariant variety (|20p tells us where the map of 
period 3 should start. Every point on (|20p is a candidate of the period 3 map. The elliptic curve 
is embedded in this invariant variety as a set of 3 points, if the initial point is on it. We can 
view this variety as a subspace of the set of all elliptic curves, which are restricted to 3 periodic 
motion. It is a highly nontrivial observation that the intersections form a surface characterized 
by certain specific relations among the invariants of the map alone. The existence of such an 
variety in any integrable map has not been known, to our knowledge, in the literature. The 
claim of our conjecture is that if there exists an invariant variety of some period, the map is 
guaranteed being integrable. Therefore the existence of the surface (|20p of period 3 is sufficient 
to guarantee the integrability of the discrete Euler top. This is true irrespective whether some 
solutions are known or not known explicitly. 

To see other conditions, let us present all expressions of (fTBj) . (fT7|) . (fT8|) after the substitution 
of©: 

71 2) = (1 + 6-6-6), 7^ = (1-6 + 6-6), 7^ = (1-6-6 + 6), (21) 

71 3) = ((1 + 6 + 6 + 6) 2 - 4(i + 66 + 66 + 66)) 

x ((1 + 6 - 6 - 6) 2 - 4(6 2 - 66 + 66 - 66)) , 
7 f = ((1 + 6 + 6 + 6) 2 - 4(i + 66 + 66 + 66)) (22) 

x ((1 - 6 + 6 - 6) 2 - 4(6 2 - 66 - 66 + 66)) , 
73° = ((1 + 6 + 6 + 6) 2 - 4(i + 66 + 66 + 66)) 

x ((1 - 6 - 6 + 6) 2 - 4(6 2 + 66 - 66 - 66)) , 
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7r 4) = ((i - 6) 2 - (6 - 6) 2 ) ((6 - i) 2 (2(1 + 6 - 6 - 6) 2 - (6 - i) 2 ) 

+ (6 - 6) 2 ((2 + 26 - 6 - 6) 2 + 466 - 86)) , 
7 f = ((l - 6) 2 - (6 - 6) 2 ) ((6 - i) 2 (2(i - 6 + 6 - 6) 2 - (6 - i) 2 ) (23) 

+ (6 - 6) 2 ((2 + 26 - 6 - 6) 2 + 466 - 86)), 
if = ((i - 6) 2 - (6 - 6) 2 ) ((6 - i) 2 (2(1 - 6 - 6 + 6) 2 - (6 - i) 2 ) 

+ (6 - 6) 2 ((2 + 26 - 6 - 6) 2 + 466 - 86))- 

The conditions (|21[) for the period 2 impose 

6 = 6 = 6 = i- (24) 

/Q\ fO\ fO\ 

The second factors of 7) , 72 , 73 vanish simultaneously only when the point (xi,X2,x 3 ) is 
on the lines defined by 

{x|6 = i n 6 = 6} u {x|6 = i n 6 = 6} u {x|6 = i n 6 = 6}- (25) 

After some manipulation we find that 7^ , 72 ; 73 also vanish simultaneously iff the point is 
on these lines. Therefore every periodic point of period 2, 3 and 4 are on the lines of (|25p . if it is 
not on the invariant variety (|20p of period 3. We now notice that the points on these lines f)25f) 
vanish the denominator of the map ([3]), or equivalently the function (1 — 6 — 6 — 6) 2 — 4666- 
From this observation we are convinced that the Euler top has no periodic point of period 2 
and 4 as long as the top is not axially symmetric, whereas the periodic points of period 3 form 
the invariant variety i>w of (|20p . 



5 Axially symmetric top 

By studying the discrete Euler top we have found an invariant variety of periodic points in 
the case of period 3. If we adopt our conjecture in Section 1, this means that the system is 
integrable. We also found that there is no periodic point of period 2 and 4 if the top is not 
axially symmetric. Our method enables us to search the periodic points of larger period. Instead 
of carrying out further the cumbersome algebraic analysis, however, we conclude this paper by 
studying the cases of symmetric top. 

If the top is totally symmetric, i.e., h = h = h, the equations Q show that (x\,X2, x 3 ) 
remain constants. This is a top which never changes its angular velocity in all directions. When 
the top is axially symmetric, such as h = h, the motion is governed by ((4j) . As we repeat the 
map n times we get 

x[ n) = Xl , 

X^ = %2 cos(nfi) + X3 sin(nJl), (26) 
X^ = X3 cos(nQ,) — X2 sin(nfJ). 

The periodicity condition of period n in this map can be read off directly from (I26p as 
cos(nO) = 1, or VL = This condition fixes the values of x\ for each period, according to the 
rule 



, 2/ 2 / l-coB(27r/n) 1 . _ 

x 1 = ±n n - -, fi n = A — — - - n=l,2,3,.... (27) 

h-h V 1 + cos(27r/n) 
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For small n's we have 

Ml = 0, ^3 = V3, [i4 = 1, ^5 = \/s - 2^/5, [i 6 = 2-V3, .... 

The conditions determine planes which are orthogonal to the axis of symmetry and intersect 
the axis at certain points defined by (|27h . different for each period. These planes 



«2L symm = {*K = ^ (j^j^ } > " = 2, 3, 4, . . . (28) 
are the invariant varieties of periodic points characterized by the relations among the invariants: 
I2H1 — H2 



1 + ^^2-/1' 



In terms of the variables (6,6,6), the conditions I2 = I3 and (|27p are equivalent to 
(6,6,6) = (-^n,°, 0)- We notice that vf] ial symm is a special case of -y (3) in d2DD - When 
n = 4 there is no invariant variety for generic values of (Ii, 12,13), hence we are not able to de- 
rive f^jaj symm as a special case. We notice that the periodicity conditions (6, 6, 6) = ( — 1, 0, 0) 

in the case of period 4 are compatible with the conditions 7^ = 73 = of (I23p . 

The meaning of the planes presented in (f28|) is quite interesting. Because of the symmetry the 
angular velocity along the symmetry axis is constant as it is expected naturally. An interesting 
feature is that the value of this angular velocity x\ is 'quantized' to some specific values deter- 
mined by the shape of the top and different for each period, such that the 'angular velocity' £1 
is quantized to 2ir/n to generate the periodic maps. This is true irrespective to the values of 

(n) 

other angular velocities X2 and X3. The generation of the invariant varieties v^ al symm follows 
to this fact. 
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